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1. Introduction 

In this work we consider supersymmetric Yang-Mills quantum mechan- 
ics (SYMQM) [1, 2]. Theory for SU(N -> oo) gauge group in D = 10 
dimensions is related to M— theory and allows researches on the BFSS hy- 
pothesis [3]. For a smaller number of dimensions and a smaller number of 
colours, like in our case where D = 4 and N = 2 respectively, such me- 
chanics pose excellent theoretical laboratory [4] to test amazing properties 
of supersymmetry, i.e. coexistence of discrete and continuous spectrum [5], 
action of supersymmetric generators or unique features of SUSY vacuum 
states [2, 6, 7, 8]. Moreover SYMQM without fermions describe glueballs 
which are also considered in many non-supersymmetric theories [9, 10, 4]. 
Here, in order to calculate the energy spectrum of our model the method 
proposed by van Baal in Ref. [8] is used. In this work eigenvalues as well as 
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eigenfunctions of Hamiltonian of the model are computed. Finally, to dis- 
tinguish the localized states form the non-localized ones the virial theorem 
is applied. 



2. Hamiltonian 

The zero- volume Hamiltonian [9, 8], i.e. in the long wave approximation 
[2, 6], reads 

H = -\Y.(-^) 2 +^E(Af ~i E ^VA^, (2.1) 

i,a 1 i,a i,a,b,d 
y v ' " „ ' v „ ' 



where in the bosonic part, i.e. in Hb, Ht is the kinetic Hamiltonian, Hy 
with ^ 

Bi = ~2 E £ ijk £ abdC b jcf, (2.2) 
i,j,k,a,b,d 

corresponds to a bosonic potential while the fermionic part is denoted by 
Hp, where = r J are Pauli matrices. Bosonic variables, cf, have colour 
indices a = 1,2,3 and spatial indices % = 1,2,3. Except the colour indices 
the anticommuting Weyl spinors, A", have spinor indices a = 1,2. Thus, in 
this system the maximum number of fermions is 3 x 2 = 6. 



3. The cut Fock space 

The operators of the number of fermionic quanta, hf, and the total angu- 
lar momentum of the system, j commute with the Hamiltonian. Therefore, 
solving the eigenproblem of the Hamiltonian we can separately consider the 
sectors with fixed values of rip and j. The system has particle-hole symme- 
try so it is enough to consider sectors with uf = 0,1,2,3. Unfortunately, 
the operator of number of bosonic quanta, ns, does not commute with H . 
Since the model has an infinite number of bosons solving it numerically we 
have to cut somewhere the Fock space off. A good choice for this cut-off is 
a maximal number of bosonic quanta in the system, B >ub- 

We are especially interested in the sector where uf = 2 and j = 0. 
In this sector the localized and non-localized states coexist. Moreover, it 
contains the supersymmetric vacuum state. To find the energy spectrum we 
use van Baal approach [11, 8]. For j = and tif = 2 the Hamiltonian can 
be rewritten in terms of three bosonic variables 



r 2 =^(c|) 2 , u = r-*Y,(Bj) 2 > v = r~ 3 detc, 

j,a j,a 



(3.1) 
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or equivalently in (xi, X2, X3) where 

r 2 = Y,*h u = r- 4 J2 x hh v = r- 3 H Xj . (3.2) 

j i>j j 

For example the bosonic potential has a form 

H v = u r 4 /2 = {x\x\ + x\x\ + x\x\) /2 . (3.3) 

The minimum of Hy is localized in six valleys along the The other 

parts of the Hamiltonian have much more complicated structure [8, 12]. 
We construct the Fock space acting and A^ on the empty state [8, 4] : 

H= E ^..c-At +1 ...Al|0>, (3.4) 

contractions 
{oi, . . . , a r } 

where sum goes over gauge and rotation invariant combinations of the op- 
erators. For uf = 2 we have two independent ways of the fermionic action: 

X/ = -2i £ £c)b , a A^(^°)^A 6 >) , J a6 = - AJAj^lO) , 

c,b,a,$ c,b,a,$ 

(3.5) 

where cr j0 = ^Tj and e aj g = e A/ g = — IT2 lowers spinor indices. Making 
contractions of (3.5) to bosons we obtain six independent invariants 

\ei(u,v)) = Ej^/^a, |e 4 («,t>)> = E^Jab, 

|e 2 (u,«)> = Ej^Bf/r^a, \e 5 (u,v)) = £ a , b c«c*/f 2 , 

\e 3 (u,v)) = Ej^Act/r^a, \ee(u,v)) = E a ,b J<* ■ 

(3.6) 

Other basis vector can be obtained acting with invariant combinations 
of bosonic variables, i.e. (r,u,v), on these six vectors. This gives following 
basis vectors 

6 

\n) = ^2h^(r,u,v)\e m (u,v)) , (3.7) 

m=l 

where h^(r,u,v) are arbitrary functions. Following Refs. [8, 12] they are 
chosen as eigenfunctions of the harmonic oscillator. In this basis the matrix 
of the kinetic Hamiltonian, Ht, is tridiagonal. 
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4. Diagonalization of the Hamiltonian matrix 

The Hamiltonian (2.1) in the basis (3.7) gives the eigenequation 

n 

where H nn = (n'\H\n) is the Hamiltonian matrix [8, 12] while are its 
eigenvalues. The matrix for cut-off B < 11 is plotted in Fig. 1. Intensity 
of gray scale shows the amplitude of the matrix elements. We can see five 
broad branches of the bosonic Hamiltonian, which corresponds to change of 
n-B by 0, ±2, ±4 and two other narrower branches of the fermionic part, Hp, 
where this change equals ±1. For B = 60 this matrix has 10416 x 10416 
elements and its diagonalization is much more time consuming 1 . 



B 34567 8 9 10 11 




Figure 1. The 112 x 112 Hamiltonian matrix in the base (3.7) for cut-off B < 11. 



The dependence of the energy, E^, on the cut-off is shown in Fig. 2. 
We can see two different behaviours of the energy levels. The first levels 
converge rapidly. They correspond to localized states. The other ones never 
converge but they slowly fall down as E(B) ~ 1/2?. Cut-off analyses of free 
particle systems gives conjecture that the latter levels gives the continuous 
spectrum at B — > oo [13]. 

1 The evaluation time grows exponentially with B. To speed up calculation the van 
Baal's Mathematica code [8] was rewritten to CH — h 
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10 20 30 40 50 60 B 



Figure 2. Energy spectrum as a function of cut-off B > ub 
The eigenfunctions of H read 

6 

\* k (r,u,v)) = 2>fcl"> = J>2 ^ C(w)|e m ) . (4.2) 

ra n m=l 

Examples of such wavefunctions were presented in [12]. Indeed, the rapidly 
converging states are localized in the centre of the x— coordinate system. 
On the other hand, the non-localized ones penetrate the system along the 
potential valleys [12]. For these states realization of the B — > oo limit is 
more complicated and it is performed with fixed energy E and changing k, 
i.e. 

\& E \r,u,v)) = lim \$ k (E K =E,B)(r,u,v)). (4.3) 

Given the explicit form of the energy eigenstates using (4.2) we can 
calculate averages of other operators 

(0} k = (* k (r,u,v))\0\<!> k (r,u,v)) . (4.4) 

The only technical problem is to rewrite the operators in terms of (r,u,v). 
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5. Virial theorem 

In order to distinguish the localized states from the non-localized ones 
we apply the virial theorem. One can derive it from the Heisenberg equation 



dF _ OF 1 . ~ 



(5.1) 



For a motion in a compact space the virial, x-p, is a limited physical variable. 
Therefore, its average does not change with time: 



= ^(x-p) = U[x-p,H}), 
at in 

where the second equality follows from (5.1). 



(5.2) 



EfHtfl4 ;i m;;!i :; ;: ; : t 



* ** 




I -3 



25 30 35 40 45 50 55 60 

B 

Figure 3. The test function (5.6) for the bound states as a function of cut-off 
B > riB for the the lowest fifteen energies. Here to guide the eye joined triangles 
correspond to states from discrete spectrum while stars are related to states with 
continuous spectrum 

In the case of Hamiltonian H = T + V where T is the kinetic energy and 
the potential energy scales as V(ax) = a n V(x) one obtains from (5.2): 



-2(T)+n(V) = 0. 



In our case 



x ■ p = —ir 



dr ' 



and the consecutive parts of the Hamiltonian scale as 

HT(ar) = a~ 2 HT{r) Hv(ar) = ce 4 Hv(r) and Hp{ar) 



(5.3) 
(5.4) 



aHp(r) . 
(5.5) 
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Therefore for the bound states the test function 



/ = -2(H T ) + A(H V ) + (H F ) = . 



(5.6) 



Applying (4.4) the averages (Hx) are calculated over the eigenstates of 
the Hamiltonian (4.2) as 



Substituting the obtained values of the averages to (5.6) we can check 
the relation (5.6) for the bound states. We plot the test function (5.6) as 
a function of the cut-off f(B) in Fig. 3. One can easily see that for the 
localized states the relation (5.6) is fulfilled with a very good approximation 
even for not so high B. 

For the states with the continuous spectrum, condition (5.6) is not sat- 
isfied. This is seen in Fig. 3 where we compare the test functions for the 
bound states and the states related to the continuous spectrum as a function 
of cut-off B > us- The test function (5.6) related to the bound states are 
nearly equal to zero while the test function for the states which form the 
continuous spectrum seem to grow with B. 



In this work we discuss supersymmetric Yang-Mills quantum mechanics 
(SYMQM) [1, 2] in four dimensions for SU(2) gauge group. Investigation of 
such models allow to understand complicated properties of the supersym- 
metric theories, i.e. coexistence of localized and non-localized states, non- 
triviality of the supersymmetric vacuum, which are common for different 
supersymmetric theories. 

We focus on the sector with the number of fermionic quanta rip = 2 and 
the total angular momentum j = 0. This sector possess both discrete and 
continuous spectrum [5], and the supersymmetric vacuum state [8, 4, 12]. In 
order to find the energy spectrum we use a method proposed by van Baal in 
Ref. [8] where the cut-off of our Fock space, i.e. B, is defined as a maximal 
number of bosonic quanta, n#. To confirm localization of states the virial 
theorem is used (5.2). 

The result of the Hamiltonian spectrum agree with the previous works 
[4, 8, 14]. However, here the calculations have been performed for the very 
high cut-off B = 60. For this cut-off not only the spectra of bound states 
converge but also the corresponding eigenstates [12]. 

Our calculation shows that the quantum virial theorem is applicable for 
the systems with more complicated potential, e.g. one considered in this 
work which consists of two different parts. Moreover, the virial theorem can 



Ex = {H x )k = {<S> k (r,u,v)\H x \<S>k{r,u,v)) . 



(5.7) 



6. Summary 
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be used to determinate localization of the states where it is not possible to 
calculate the Hamiltonian eigenstates and where the averages of operators 
can be computed in a different way. Thus in some cases the virial method 
as the localization test of states is easier and more applicable that direct 
calculation. 

In the future using this method one can calculate matrix representations 
for other operators (4.4) and test various properties and laws for this model. 
Furthermore, using similar methods one can also try to solve the models in 
more dimensions and for different gauge groups [15]. 
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